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INTRODUCTION 


POSSl  is  a FORTRAN  code  for  computing  the  response  of  an  axl- 
symmetric  soil/structure  interaction  problem  in  cylindrical  (r,  z) 
coordinates.  The  acronym  POSSl  stands  for  Problem  of  Soil  Structure 
Interaction.  The  basic  two  dimensional  configuration  of  the  prol)lem, 
consisting  of  a lined  cylindrical  cavity  of  finite  length  extending  from 
the  surface  into  a semi-infinite  half-space  and  loaded  Internally,  is 
shown  in  Fig.  1.  Modification  of  certain  parameters  in  the  program  do, 
however,  permit  use  of  the  code  for  one  dimensional,  or  free  field 
calculations. 


Development  of  the  computational  approach  for  POSSl  has  been  pre- 
sented in  Refs.  [l]-[4].  This  report  is  the  final  one  in  that 
series,  and  presents  not  only  information  about  use  of  the  code,  but 
also  adds  certain  theoretical  features. 

The  structure  which  lines  the  cavity  is  considered  as  an  axlsymmetrlc 
thin  cylindrical  shell  with  bending  stiffness,  welded  to  an  axisymmetric 
circular  bottom  plate.  In  the  previous  reports  the  plate  was  considered 
to  have  bending  stiffness  only.  The  work  presented  in  this  report  Includes 
extensional  motion  of  the  plate.  Also  presented  is  the  development  of 
equations  governing  the  maintenance  of  a right  angle  at  the  joining  of 
plate  and  sliell,  see  Fig.  1.  In  addition  to  the  equations  for  these  new 
features  the  basic  equations  given  in  Refs-  ( M.  I I liavi-  been  repeated  here 
for  convenience,  and  grouped  according  to  the  order  of  computation  in  the 
POSSl  code. 
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The  material  surrounding  the  structure,  and  interac-ting  with  it. 


is  considered  to  be  a bilinear  material,  with  hysteresis  in  tlie  response 
of  volumetric  effects,  but  linear  response  in  shear.  Provision  is  made 
for  horizontal  layering  of  the  material,  with  different  material  properties 
for  each  layer. 

The  computational  scheme  used  to  describe  the  response  of  points  in 
the  field  surrounding  the  structure  is  the  pseudo-characteristic  scheme 
combined  with  a fractional  step  method  that  was  developed  in  Refs.  [ll-[^l. 
As  in  Ref.  [3]  this  scheme  will  be  hereinafter  referred  to  as  the  PC  scheme-. 

The  governing  equations  for  the  motion  of  the  structural  lining 
were  obtained  from  a method  of  the  variational  calculus  applied  to  the 
finite  difference  form  of  the  total  strain  energy.  The  right-angle  condi- 
tion at  the  joining  was  Introduced  as  an  equation  of  constraint  with  the 
addition  of  a Lagrangian  multiplier.  The  resulting  equations  are  first 
order  finite-difference  in  nature. 

A number  of  options  are  built  into  the  POSSl  code.  They  include: 
location  and  type  of  applied  load;  dimensional  or  non-dimensional  output; 
lined  or  unlined  cavity;  right-angle  at  plate-shell  connection  maintained 
or  not  maintained;  horizontal  layering  in  the  field,  or  no  horizontal 
layering;  graphic  or  printed  output,  or  both;  one-dimensional  configurations 
in  either  r-t  or  z-t  coordinates  as  special  cases;  the  shell  and  tlie  plate 
may  have  different  material  properties  or  thicknesses.  These,  and  otlier 
features,  are  discussed  in  Section  111  of  this  report. 
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In  ordor  to  illustrate  use  t)l  tlu‘  I’OSSI  rode  typical  output  from 
three  sets  of  eompuLat iiuis  are  presented.  The  first  ease  is  lOr  an 
unlined  cavity  in  a layered  material,  subjected  to  a time  deiu'iuient  , 
and  location  dependent  internal  pressure  ii(z,  t).  Two  one-dimensional 
subcases  are  also  very  briefly  discussed.  Tlie  second  c.ise  is  for  .1 
structurally  lined  cavity,  in  a layered  material  with  a right  any.le 
maintained  at  tlie  sbell-|)late  connection,  subjected  to  an  intern.il 
pressure  p(z,  l).  Results  in  this  case  are  presented  in  non-dimensional 
form.  The  third  case  is  a rerun  cif  the'  sc'cond,  witli  results  prc'scnited 


in  dimensional  form. 


A . Field 


The  material  surrounding  tlie  structure  acts  linearly  in 

shear  (i.e.  C,  the  shear  modulus,  remains  constant),  but  bilinearly 

in  volumetric  response.  That  is,  the  bulk  modulus  K has  two  possible 

values:  p virgin  loading,  K = for  (-.1^)  less 

than  or  equal  to  the  previous  compressive  maximum.  Figure  2 shows  how 

this  restricts  the  dependency  of  > on  .1 , . The  model  continues  onto 

V 1 

the  tensile  side  of  the  Jj^,  curve  along  a straiglit  line  of  essentially 
unloading  slope,  unless  tension  is  experienced  initially,  in  which  case 
the  loading  slope  is  used.  .See  Fig.  2. 


The  basic  differential  equations  which  describe  the  material  motion 
are  the  same  as  those  for  a linear  elastic  continuum.  Only  the  parameters 
K (bulk  modulus)  and  v (Poisson  ratio)  must  be  appropriately  chosen, 
according  to  the  criteria  mentioned  above,  for  conditions  of  loading  or 
unloading.  G (the  shear  modulus)  remains  constant  throughout  the 

computation.  In  non-dimensional  terms  these  equations  are 

9V  9o  a - a„  Set 
r _ r r 0 rz 

9x  9r  r 9z 

9V  9a  a 9o 
z rz  rz  ^ 

^ 9x  " 9r  r 9z 

90  „ 9V  V 9V 

r \2rr.  V ,r.  z,, 

97  ^ 9r  1 - V ^'r  + ^ 1 

9a  „ 9V  V 9V 

-9T=''^  t-9f  <7-^-9T>) 


9n  - 9V  9V 

= ma2  f-^  + --'■-I 
dr  s dr  d z 
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Till'  symbols  aro  ileliiu'd  ia  Table  1,  as  are  all  I be  parameters, 
tile  stresses,  vt'loeilies  aiul  d is))  I aei'iiietU  s wliieli  a|)pear  in  t lu'  basii- 
eiiiiaLions  for  both  field  and  struetiire.  Tbe  (|iiantities  appearing  in 
fable  1 are  given  in  both  d imi'iis  i ona  I and  non-dimensional  form. 

Using  the  method  of  fraitional  steps  the  basic  equations,  with  r, 
z anil  t dependi'iic  i es  are  separated  into  two  problems:  in  r,  t coordinates; 

in  /,  t coord  i nat  I's . Kacb  is  tlu’n  a oiu'-d  imens  iona  1 wave  propagation 
problem  with  its  own  si't  ol  1’-,  S-  and  "zi'ro"-charac  t er  i s t i c s . Integration 
is  set  lip  along  tlii'se  characteristics  resulting  in  the  equations  which 
are  summarized  in  fable  2.  (file  superscripts  refer  to  points  in  the  spa.  t 
time  meshes,  shown  in  I'igs.  la  and  lb.) 

In  order  to  obtain  values  of  stress  and  velocity  at  points  intiT- 
mediate  to  the  mesh  points  (i.e.  at  the  R and  Q pts.  shown  in  Tigs,  la 
and  lb,  between  I.-l,  1.  and  L,  l.+  l,  respectively)  a three  point  inter- 
polation scheme  was  used  leading  to  a "second  order"  inti'gration  scheme 

2 2, 

which  I'xhibits  errors  of  Ar^  ami  Az  in  r-  and  z-directions  respect  i ve  1 v . 

As  an  example  of  this  interpolation  scheme,  if  F , is  the  value  it  f 

'l> 

the  quantity  to  be  found  from  known  quantities  F^  , ''|+] 

t hen 

(.-1)1^^,  I ' mm, -I  (..) 
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TABLE  I PARAMETERS  AND  NON-DIMENSIONAL  PARAMETERS 


i 

PARAMETERS 

NON-DIMENS lONAL  PARAMETERS 

j 

COORD  1 NATES 

1 

i ^0 

Reference  radius.  Usually  the  nominal 
radius  of  the  cavity 

r = 

~ radial  coordinate^ 

^ C’  ” cavity  radius 

i 

R. 

i 

i 

Z Radial  and  depth  coordinates 

z = 

Z ! 

depth  coordinate  ' 

“ i 

1 

1 t 

Time 

i 

1 ^ 

Reference  wave  speed 

i = 

ct  . ! 

~ time  1 

*^0 

i" 

1 

Length  of  the  cavity 

1 = 

cavity  length 

0 

PARTICLE  VELOCITIES  AND  DISPLACEMENTS 

1 

ijcii 

iL 

1 

Radial  particle  velocity 

V 

r 

ecu 

u 

z 

Longitudinal  particle  velocity 

V 

Z 

Po’ 

pc^u^ 

u 

r 

Radial  displacement 

W 

r 

■ 'o'’o’  ! 

1 

OC^U  ^ 

u 

Longitudinal  displacement 

W 

z 

= wf  ! 

1 

f AND  STRESSES 

I Pg  Reference  stress  from  applied  load 

o Radial  stress 
r 

<1  Lonj;i  t ud  inal  stress 

7. 

Hoop  stress 
I ’ Slu-ar  stress 

X7. 


i 


I 
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PARAMETERS 


NON-DIHENSIONAI,  PARAMETERS 


First  invariant  of  stress 
= "r  ^ "z  ^ % 


jN  , N Eongitudinal  forces  in  the  shell 
j ^ P 

and  plate,  respectively 


0,0  Shear  forces  in  the  shell  and 

I s p 

! plate,  respectively 


1 Pr 


N 

N 

s 

M — 

_.2. 

t 2 ^ 

N — 

p 

^ 2 

pc 

pC 

Q 

Q 

s 

. -P  . 

% ■ 

,10“/" 
''  7 

pc 

0C“ 

|M  , M Bending  moments  in  the  shell  and 

I s’  p 

I plate,  respectively 


f)C“ 


, M^  = — i>-—  1 


i ^lATER  l AF,  P_RO?.!^^*lTUc’ 
p Reference  density 


(1 . Material  density  of  field  (a  different 

p.  will  be  specified  for  each  layer) 


^i 


C Reference  wave  speed 
K + - r 

'2  LD  3 , 

C.  rx  = liOacling  P-wavo  spoea 

i laJ  P 


.2  %N  3 ^ 


Unloading-Reload ing 
P-wave  speed 


El)  C 

V = 

UN  C 


2 , 

C = Shear  wave  speed 
s p 


'"'.'I 

- 

20 

'’^ED 

"+ 

Yc 

’"t'N 

- 

2C 

V 

Yc 

Eoading  Poisson  ratio 


Poisson  ratio 


' ^ 1 D’  ^UN’  ^ specified  for  eacli 


I 


PARAMETERS 


non-dimf:nstonal  parameters 


. STRUCTURE  PROPERTIES 


y2  ^ ^ 


E 

Op  i 


Fundamental  frequency  for  shell  and 
plate,  respectively 


li  , h Shell,  plate  thickness 
s p 

I h 2 , h 2 

“» ■ n <-sJ>  ■ “p  ■ T2  <i^> 


Q R„ 
= s 0 

s C 


^ R„ 

‘>n  = 

P c 


h h 

R„  ’ R 
0 0 


a , a 
s p 


B = — 6 = ■ - 

s p h ’ p p h 
s s p p 


3 , B 

s p 


Eh  Eh 

D = , 0 = 
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TABLE  2 EQUATIONS  FOR  COMPUTATION  OF  STRESSES  AND  VELOCITIES  IN  THE  FIELD 


CHARACTERISTICS  FOR  THE  r- DIRECTION 
FIG.  3a 


XXX 


FIG.  3b 
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where  the  appropriate  value  of  A determines  nt  as 


In  the  POSSI  Code  this  interpolation  scheme  has  been  used  at  all 
points  in  the  field  in  both  the  r-  and  z-directions. 


The  PC-schcme  as  outlined  above  does  not  provide  for  following 
or  maintaining  shock  fronts.  (Computational  results  from  a loading  with  a 
very  short  rise  time  would  therefore  be  dispersed. 


B.  Structure 


The  structure  lining  the  cavity  Is  a cylindrical  shell  with 
bending  stiffness,  welded  at  tlic  bottom  to  a flat  circular  plate  which 
deforms  both  by  bending  and  extension.  Tlie  equations  of  motion  in  both 
radial  and  longitudinal  directions  for  botli  structural  components  have 
been  obtained  by  a variational  procedure  applied  to  the  total  strain 
energy  expressed  in  finite  difference  form.  Tlie  expressions  for  potential 
energy  of  tlie  shell  (including  bending  energy,  extensional  energy  and 
work  done  by  the  loads)  plus  the  kinetic  energy  of  the  shell  are  given  in 
Appendix  B of  Ref.  [A ] . Repeated  liere,  for  convenience,  in  terms  of  the 
non-dimensional  quantities  defined  in  Table  1,  the  total  energy  for  the 
shell  is 


trD  R„  R,,  , ,,  DW  ^ „ ;iW 

If  = ( .9  ) { 1 [(  /:)  + + 2v  W — ''l 
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0 a 2 T 
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A similar  oxpross  i on  Tor  t lii’  plato 
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iiu'liuli's  t ho  oiuTgy  ol  oxli'iis  iona  1 plato  mot  ion  which  was  oxoliuloil  I rom 
the  plato  oiuTgv  t'xprossion  given  in  Ri’f.  | •♦  | . The  ilorivativos  in  I hosi’ 
oxprossLt)iis  arc  written  in  f i n i t o-d  i I I oriMU'o  li'rm  at  oviMv  point  on  t lie 
shell.  Central  ilin'erenees  .iri'  nsi-d  lor  non-ln'iinda  rv  points,  rorward  or 
backward  d i I' fe  ri-nees  , as  approi>r  i at  i\  a re  used  tor  points  in  the  no  i I'.lilu' r- 
hootl  of  the  corner.  A v.i  r i.i  t i ona  I proeednri'  is  then  applied  to  the  I ot  .]  I 


I'nergy,  first 

with  respect  to  r, 

ml  i ;i  1 tl  i sp  1 acement  , 

t hen  with  re 

sfiec  t t i' 

1 ong i t ml  i na 1 

d i sp  1 aci'ment  . See 

Ainn'mlix  It,  Ref.  | 

■t  j . 'fhe  resn 

It  i t wo 

etpi.t  t ions  o 1 

mot  ion , one  i n t he 

r-tl  i ret' t i tni , one 

i n t he  /.-tl  ire 

c t i t'n  at 
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The  snhscript  "m"  is  writti'n  as  r , , r , depeiulini',  upon  the  dire.  l i 

s s p p 

ol  mot  ion  lu’ing  computed,  and  whether  t lu'  term  ci'ini's  I fi'.ii  shell  oi  plali' 

elli'cts.  Tlu'  (lu.intities  a'  art' 

m 

h I,  2 

A = f,  o + .1  IP.’ZWZ  + \)  1)1  ZWK  + , I'tZ.KKi 

z s rz  s s t; 
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I,  I L ^ I 
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tor  llu*  sliell,  .liicl 


for  the  plate.  The  siiiierscr  ipt  "I,"  refers  to  a point  on  the  striieture, 
see  Fig.  4.  Tiie  superscript  "N"  refers  to  tlie  new  value  of  the  quantity 
at  point  L.  The  symbol  D2RWZ  stnmis  for  the  second  derivative  with 
respect  to  r of  W^,  IllZWR  stands  for  tlie  first  derivative  with  respect 
to  7.  of  W^,  etc.  The  finite  difference  expressions  for  these  symbols 
are  tabulated  in  Table  3a  and  3b  (for  the  shell).  Tables  4a  and  4b  (for 
the  plate). 

Kquations  (10)-(14)  are  expressions  for  velocities.  The  displace- 
ments are  found  as 
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'the  POSSI  code  also  provides  for  computation  of  forces  and  moments 
within  tlie  structure.  For  the  shell  portion  the  normal,  shear  and  moment 
(]uantities  per  unit  of  circumferential  length  are 
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NUMBERING  OF  POINTS  ON  STRUCTURE 


TABLE  3a  E'^LATIONS  OF  MOTION  FOR  THE  SHELL 


Kq . i'll)  for  definition  of 


TABLE  3b  EQUATIONS  OF  MOTION  FOR  THE  SHELL 
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TA3LL;  4a  EQUATIONS  OF  MOTION  FOR  THE  PLATE 


*See  Fi;^.  4 

**  See  Eq.  (27)  for  definition  of 


I'or  the  plate  portion  the  normal,  shear  and  moment  quantities  per  unit  of 
e.ircumferential  length  are 
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See  Table  1,  item  4 for  definition  of  the  corresponding  non- 
dimensional  and  dimensional  quantities. 
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1.  Shf  1 I - 1’ 1 a t c*  (aMiiu'c  t i (,)ii 


H(|iiatioiis  (10)-(14)  dose  r i bo  volocities  at  structure  points  provldiuK 
no  restriction  is  made  concerning  the  relal  ive  angli'  between  tlie  shell 
wall  andt'ac  plate  bottom.  It’  a right  angle  is  to  be  maintained  at  this 
joining  then  the  equation  is  modified  for  points  near  the  corner. 


The  I'ondition  of  a right  angle  at  the  shell-plate  interface  is 
stated  by  the  constraint  condition 

■)W  3W 


the  subscript  "LC"  defines  the  point  at  the  interface,  see  I'ig.  4.  This 
is  introduced  into  the  variational  procedure  by  the  Lngrangian  multiiilier 
A such  that  the  total  variational  expressions  are 


d(r  + u ) 


where  W indicates  displacement  at  point  "l."  i n e i t her  the  r-d  i rt'c  t i on  . 
"'h 

or  the  z-direction,  and  U , I'  are  the  total  energv  expressions  given 

s p 

in  K(|s.  (8)  and  (9).  For  convenience  the  energy  variat  ion  is  wr  i t t ini 
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[where  t hi'  i|uantities  are  given  in  K(|s.  ( 1 1 ) - ( I 4 i | and  l lie  constraint, 

Kq  . (17),  i writli'n  in  f i n i 1 1'  difterence  form 
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Thc'n  Eqs.  (18),  (19)  and  (20)  combine  to  give 
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A fifth  equation,  which  allows  for  solution  of  A is  obtained  from  the 
second  derivative  with  respect  to  time  of  Eq . (19) 
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Combining  Eqs.  (20)-(24)  yield 


A = 


. ,.LC  ,LC  ,LC-1,  , ,LC  IX  ,LC  JX+1 

Ar(A  + A - A ) + A + A + A - A 
c*  r r r z z z z 

Az s .P  . s s p p p_  _ 

' ' 9 ' T 


.Ar,“ 


Y-x-  fi  + (r.)  1 + (-.-)  + r 


1 + 


V\z 


',Z 


1 


Y r 


IX_ 

IX+1 


where 


9 

JI 

■ z 

•,) 

P 


2 '"iX 


See  Fig.  4 for  points  corresponding  to  superscripts  L(\  l.t'H,  I.C-I,  etc 
In  a layefi’d  medium  it  is  understoo<l  tli.it  t hi'  'z  term  .ippe.iring  in 
eqii.it  ions  ri'l.ited  to  the  m.i  i nt  en.ince  of  .i  right  , ingle  at  the  corner 
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(21) 

(22) 

(23) 

(24) 

(25) 
(2b) 
(27) 


slumld  be  the  value  of  the  Az  in  that  layer  in  wliicli  tlii'  plate  sits. 


■I 


For  a structure  in  which  there  is  no  right  angle  restriction  at 
tht'  corner  the  equations  of  motion  for  all  points  are  given  bv  Fq . (10). 

For  a structure  with  a right  angle  rest  rii  t ion  Fc|s.  (21)-(24)  must  he 
used  as  equations  of  motion  lor  points  1,C-1,  L(-'  and  I.C+l.  All  othi-r 
points  are  described  as  before. 
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Ill  THK  POSSI  CODE 
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A.  Conditions  for  use 

Amons  the  options  built  into  the  Possi  Code  is  a choice  of  location 
and  type  of  loading.  There  are  three  possibilities;  loading  on  the  upper 
boundary,  see  Fig.  5a;  loading  on  the  interior  of  tlie  cavity  (or  structure 
if  it  is  present).  Fig.  5b;  loading  on  both  the  upper  boundary  and  tlie 
cavity  interior.  Fig.  5c.  In  each  case  a time  variation  may  he  combined 
witli  an  appropriate  space  variation,  as  shown.  The  input  parameter 
controlling  which  of  the  tliree  loading  types  is  to  be  used  is  TLOAD.  See 
card  #2,  Table  7.  The  actual  loading  function  is  to  be  input  directly 
into  program  functions  APSGZ  and  SIGAP  as  described  in  the  next  section, 
if  no  load  is  applied  on  eitlier  the  upper  or  cavity  boundary  conditions 
of  zero  stress  are  assumed.  It  should  be  noted  that  boundary  conditions 
at  maximum  grid  lines  are  the  type  of  transmitting  boundary  obtained 
from  mirror  image  values  of  stress  and  velocity.  Gonditions  along  a line 
through  the  origin,  but  below  the  cavity  are,  of  course,  those  of  zero 
velocity  in  tiie  radial  diret-tion. 


The  other  options,  such  as  non-dimensional  versus  dimensional 
computations,  number  of  layers,  lined  or  unlined  cavity,  right  angle  at 
shell-plate  interface,  graplilc  or  printed  output,  are  also  controlled 
by  input  parameters  which  are  listed  in  detail  in  Table  7,  Appendix  A. 


The  options  of  one  dimensional  (r-t  or  z-t)  computations  are 
activated  by  choice  of  grid  dimensions.  For  a one  dimensional  r-t  compu- 
tation  choose  a lotig  slender  cavity  with  loading  on  the  interior  vertical 
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FIG.5a  LOADING  ON  UPPER  BOUNDARY 

, I 
i { 


FIG.  5b  LOADING  ON  INTERIOR  OF  CAVITY 


1 I 1 I i i I I I I 


FIG.  5c  COMBINED  LOADING 


LOADING  OPTIONS  FOR  POSSI  CODE 
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ehoiise  ;i  p(t)  lo.icli.ng  on  the  upper  boundary,  a long  slender  cavity  with 
no  loading  and  a deep  grid. 

The  numbering  system  used  to  identify  points  in  tiie  field  and 
points  on  the  shell  is  sliown  in  Fig.  b.  Sjiecial  note  sliould  be  m.ide  of 
tlie  fact  that  field  points  are  numbered  (in  the  r-direc't  ion)  from  the 
first  point  away  from  tlie  cavity,  rather  than  from  the  origin,  except  at 
points  below  the  cavity  where  numbering  dcH*s  begin  from  tiie  origin. 
Regardless  of  grid  dimensions  chosen  there  should  be  at  least  3 points 
on  the  vertical  boundary  of  tlie  cavity  and  at  least  5 points  on  tlie 
horizontal  boundary  from  r=0  to  ^ at  tlie  bottom  of  tlie  cavity. 

Without  precise  stability  criteria,  it  is  recommended  that  time  and 

X At 

space  increments  be  chosen  such  that  the  quantity  well 

below  1.0  when  there  is  no  structure,  and  below  when  a shell-plate 
structure  is  present. 

Typical  compile  time  on  a C[)C  6600  is  approximately  11  seconds. 
Typical  run  time  for  a 60  by  40  point  grid,  with  three  horizontal  layers 
and  no  structure,  running  300  time  steps  with  printed  output  every  10th 
time  step,  and  output  for  time  history  plots  every  2nd  time  step  is 
about  380  seconds.  Typical  run  time  for  a 60  by  40  grid  with  one  layer 
and  a structural  lining,  running  200  time  steps  with  the  same  fretpiencv  of 
out  put  as  above  is  about  220  seconds.  Additional  outfuit  will,  of  course, 
increase  running  times  noticeably. 

It  may  also  be  noted  here  that  all  input  is  on  TAI’K  2,  all  printed 
output  is  on  TAPF  3,  and  all  output  for  graphing  is  on  TAPK  4. 
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TYPICAL  NUMBERING  OF  POINTS 
ON  THE  STRUCTURE 


i ^TYPICAL  NUMBERING  OF  POINTS 

' / IN  THE  FIELD 
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B.  DESCRIPTION  OF  SUBROUTINES 


The  following  describes  the  essential  function  of  eacli  subroutine 
in  the  POSSI  Code.  The  Subroutines  are  listed  in  ttie  order  in  which  they 
are  compiled  in  the  code 

RZ2D  - MAIN 

Increments  time.  Activates  integration  first  in  the  r-direction 
(CALL  RDIR),  then  in  tlie  z-directlon  (CALL  ZDIR) . If  a structure  is 
present  witli  a right  angle  restriction  at  its  corner,  RZ2D  - MAIN 
corrects  tlie  shell  corner  motion  for  that  restriction.  Integrates 
velocities  into  displacements,  Eqs.  (15),  (16)  at  grid  points  on  the 
structure.  Calls  for  check  of  loading,  unloading,  reloading  conditions 
at  all  field  points  (CALL  TEST).  Initiates  graphic  (CALL  GRAPHS)  or 
printed  (CALL  SLOT)  output  as  required. 

LOCATE 

Computes  arrays  LR,  LRMX  which  control  the  extent  of  integration 
in  the  r-directlon,  and  arrays  LZ,  LZMX  which  control  tlie  extent  of  in- 
tegration in  the  z-direction.  At  eacli  time  step  the  value  of 
2 2 

V = r ^ z computed  at  each  field  point  witti  coordinates  I,  .1. 

As  long  as  a velocity  V .001  is  encountered  the  arrays  are  set  as 
LR(1)  = J + 1,  LZ(J)  =1+1.  The  values  LRMX,  LZMX  are  then  sot  as 
LRMX  = max  LR(1),  LZMX  = max  LZ(J).  Tills  criterion  was  selected  for 
disturbance  expanding  outward  from  the  cavity,  and  for  essentially  non- 
dimensional  computations  where  V = 1.0.  For  other  shape  disturbances 
or  for  some  dimensional  computations  where  very  large  or  verv  small  Vs 
are  expected  this  criteria  should  be  changed. 
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DATA 


Requests  input  data  according  to  specifications  listed  in  Table  6. 
Initiates  computation  of  constants  (CALL  CONST).  Outputs  all  pertinent 
computational  parameters  as  shovm  in  the  sample  computations.  Section  IV. 

CONST 

Initializes  arrays.  All  stresses  and  velocities  are  set  to  zero 
Computes  constants  CNRl. . .CNR12,  arrays  CLD,  CS,  CUN  as  defined  in  the 
list  of  common  variables  in  Appendix  C.  Computes  other  constants  from 
input  data  which  are  basic  to  the  code. 

MIR 

Integrates  equations  of  motion  in  the  r-direction,  at  all  field 
points,  by  evaluating  the  upper  set  of  equations  shown  in  Table  2.  First 
integration  is  done  in  Region  1 as  shown  in  Fig.  7a.  Then  Subroutine 
SHELL  is  called  at  the  end  of  each  I*"^  line  of  integration  to  perform 
the  Integration  of  the  shell  equations  of  motion  in  both  the  r-  and  z- 
directions,  (if  there  is  a structure),  or  to  compute  the  motion  of 
field  points  at  an  applied  pressure  boundary,  if  there  is  no  structure. 

If  a structure  is  present,  all  its  mass  is  considered  to  be  distributed 
along  the  line  defining  the  cavity,  and  all  motions  of  points  on  that 
line  are  goverened  by  the  structural  equations  of  motion.  If  no  structure 
is  present  the  line  defining  tlie  cavity  is  all  material  and  requires 
integration  of  the  field  equations  in  both  the  r-  and  z-directions . 
Therefore,  in  the  case  of  no  structure,  r-direction  integration  in  Region 
1,  Fig.  7a  is  followed  by  r-direction  integration  along  the  bottom  boundary 
of  the  hole,  (CALL  PLATR) , which  is  then  followed  by  r-direction  integration 
in  Region  2,  Fig.  7a. 
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/.DIR 

Integrates  equations  of  motion  iti  the  z-direction  at  all  field  points 
by  evaluating  the  lower  set  of  equations  shown  in  Table  2.  First  the 
integration  is  clone  in  Region  1,  Fig.  7b,  then  the  integration  is  done 
in  Region  2.  At  the  t-nd  of  each  JO^*'  line  of  integration  in  Region  2 
Svihrout  ine  PLATE  is  called  to  provide  the  integration  of  the  plate 
equations  of  motion  in  both  the  r-  and  z-directions  (if  there  is  a 
structure),  or  to  compute  the  field  integration  at  points  on  an  applied 
pressure  boundary,  if  there  is  no  structure.  Then,  if  no  structure  is 
present,  for  the  reasons  stated  above,  integration  of  tlie  field  equations 
along  the  vertical  cavity  boundary  is  initiated  by  calling  svdirout  iiu>  SHKLZ. 

IIPBNI) 

Integrates  the  equations  of  motion  in  tlie  z-direction  at  field 
points  along  the  upper,  z = 0,  boundary.  Calls  fuitction  APSCZ  for 
applied  stress  boundary  condition. 

SlillLL 

Evaluates  the  eipiat  ions  of  mot  ion  for  the  sltel  1 , as  given  in 

Tables  3a  and  3b,  if  a structure  is  present.  If  tlie  cavitv  is  unlined, 

the  field  equations  of  motion  are  integrated  in  the  r-direction  lor  an 

applied  stress  boundary  condition.  Tlie  tunction  SIC.AP  is  called  tor 

application  of  o to  the  inner  surtace  ot  the  cavitv  or  slu'll,  ('  =0 

' r rz 

is  always  assumed  on  this  surlaci'). 

S112.LZ 

Integrates  the  fii’ld  e, pi. it  ions  ot  motion  in  the  ’-d  i roe  t i on  .along 
till"  vertie.il  bound. irv  of  the  c.ivitv.  This  subroutine  is  .ictiv.ited  only 
if  the  c.ivirv  is  unlined. 


Pj.A'fK 

Kv.'i  1 ii.il  I's  lIr>  ii'iis  c>l  motion  for  1 lu'  I'l.ito  as  givon  in  ral)li>s 

4a  ami  4b,  if  a st-rucluro  is  prcsoiU.  If  tlio  structuik'  is  tinliiu’d,  iho 

fiold  eqnat  ions  of  motion  are  ovaliiatod  in  t ln‘  .-t-diroit  ion  tor  an  a|iplii'd 

stress  botintlary  condition.  Tlic  function  SKIM’  is  calli'd  for  an  apt)licd 

stress  in  the  /-direction  on  tlu'  inner  surfaci’  of  the  plate  or  tlie  iinlini'd 

cavitv.  (Again  O = 0 is  assumed  on  tliis  surface.) 
r/. 

l’I.ATK 

Integrates  the  field  equations  of  motion  in  the  r-direction  .along 
the  horizontal  boundary  of  the  cavity.  This  subrout  ini'  is  .activated 
only  if  the  cavity  is  unlined. 

COKFKS 

Sets  the  coefficient  .arrays  C,  and  CS  to  the  .appropr  i a 1 1'  values  for 
each  layer.  Sets  the  coefficient  arr.ay  C .at  each  field  point  to  loading 
or  unloading  v.alues  as  determined  by  the  v lue  of  the  indicator  KOilKS . 

TKST 

At  e;ich  point  in  the  field  the  nev  v.il'ie  of  t hi-  first  invariant 
Adi  is  tested  against  the  last  previous  compressive  m.aximum  .A.llM.X.  Depending 
upon  the  outcome  of  this  test  t hi-  indicator  KODKS  is  given  a numerical 
value  to  indic.ate  conditions  of  lo.ading,  un  1 oad  ing/ re  1 o.ad  ing  .as  listi'd 
in  Table  7,  and  .‘diown  graphica  1 1 y i n F i g . 8.  Providing  for  three  cvcies 
of  load/unlo.ad  w.ls  sufficient  for  the  compiit  .at  ions  done  so  f.ar.  The  usi'f 
may  fiml,  for  certain  loadings,  further  cycling  should  be  provided  for. 

NWAl.S 

The  ev.ilu.it  ion  of  equ.at  ions  given  in  Table  2 involves  upd.it  ing  of 
qu.antities  at  o.ach  field  point.  I,,  from  non  updatcil  (pi.mtities  at  l.-l  .ami 


l.+l.  As  the  integration  proceeds  forward  in  either  tlu'  r-  or  z-dircetion 


t lio  iipdatocl  quantities  at  "L"  are  stored  in  a set  of  temporary  arrays 
(rSR,  TS/.,  rSRZ,  TVR,  TVZ,  TAJI)  until  the  updated  quantities  at  h+1 
are  evaluated  and  stored  temporarily.  Tlien  this  subroutine  NUVALS 
places  the  tem|)orary  arrays  from  point  L into  tlie  final  arrays  (Sf.R,  SC, 7., 
SGRZ,  VR,  VZ,  AJl)  of  stress  and  velocity  at  point  "L".  This  flip-flopping 
minimizes  the  amount  of  memory  required  in  the  integration. 

SllhNU 

Performs  tiie  same  function  that  NUVALS  does,  only  for  quantities 
being  integrated  along  the  cavity  boundary  in  SHELZ  and  PLATR. 

SI  GAP 

Kunction  which  sets  tlie  value  of  applied  normal  stress  on  tlie  interior 
boundary  of  the  cavity  or  the  structure,  as  a function  of  time  and  location. 
Applied  pressure  is  to  bo  considered  positive  in  this  function.  Each 
time  a different  a[>plied  pressure  function  is  to  be  considered,  the  state- 
ments in  tliis  function  must  be  changed. 

APSGZ 

Eunct  ion  witich  sets  the  value  of  applied  normal  stress  on  the  upiier 
boundary  , Z=0.  See  comments  under  SIGAP. 

SI.OT 

Directs  printed  output,  of  I'omput at  ions , samples  of  which  are  shown 
in  Appendix  H.  AM  oiiipiit  I com  Sl.Of  if;  on  iAPK  t. 

tilUPIlS 

Aiu-iimulates  output  for  later  use  hv  plotting  siibri'ut  ines . All  output 
t rom  GRAI’IIS  is  on  TAPI'  i. 
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Table  5 - Indicators  for  Loading,  Unloading  and  Reloading 


COMPUTATIONAL 

KODES  CONDITION  PARAMETERS 
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Initial  virgin  loading 

C 

Pld 

’ S ’ Pld  ' 

10 

Unload  from  first  compressive 
J,  maximum,  or  reload  up  to 

c 

Pun 

’ *^s’  ^UN  1 

^ (1) 
that  value  (--Ij^)'' 

max 

11 

Second  virgin  loading  (“J^)  ^ (“-1^/^^ 

max 

c 

Pld 

’ ^'s’  Pld 

100 

Unload  from  second  compressive 

c 

n 

’ ^s’  ^UN 

J. , or  reload  up  to  that  value 
max 

Pun 

111 

Third  virgin  loading 

c 

p 

’ ^'s’  Pf.D 

LD 

1000 

Unload  from  third  compressive 
or  reload 

C^ 

PflN 
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I 

I A.  Case  with  1 aye  red  mi'clium,  ao  s t r lu- 1 u re  . 

(Computations  were  done  for  t lie  I h rei’"  i aye  reil  nu’dium  with  no 

struetnre,  sliown  in  I'ig,  9.  The  material  parameters  for  this  ease  are 

tile  same  as  tliose  used  for  computations  in  a tliree  layeri'd  medium  with 

an  uni  Lned  cavity  di.scusst'd  in  Ref.  13).  llu'  loading  on  the  interior 

of  tlie  cavity  is  a triangular  pulse  in  t inii’ , i'ig.  10a,  hut  in  tliis  case 

-.2/. 

lias  an  exponential  variation  of  e (i'ig-  lOli)  as  opiiosed  to  the  space 

independent  pulse  tliat  was  used  in  Ref.  13 1 . Quantities  were  computed 
in  non-dimensional  form,  and  graphic  i>uiput  was  obtained  for  r.idial 
velocity  nea  r the  mi  d-po  in  t of  eticli  layer  at  r = 2.0. 

A sam|)le  of  portion  of  POSS I code  outtiut  is  shown  in  I'ig.  11.  hisled 
in  this  figure  are  all  the  input  quantities  recjuired  for  the  com|>uta  t i t'li 
of  Case  A ntsulLs. 

Radial  velocity  time  histories  at  ptiinls  msir  the  mid-point  ot  e.ich 
layer,  ;it  r = 2.0  ari-  shown  in  Fig.  12.  The  relative  li  i I fereneis  in 
arrival  times  are  of  course  dui'  to  the  ditlerent  loading  speeds  in  the 
thrt'i.'  layiTs  ( 1 . ()/ 1 . 09  1 / 3 . 1 (|  3)  . I'he  ileireasing  magnitude  ol  vi' I oc  i t \' , 
with  ik'pth  represi’uts  the  response  to  the  spatial  disttihnt  imi  ol  the 
loailing  I unction.  The  relative  length  of  t imi'  ol  a lull  pulse  in  eai  h 
laver  is  diu’  to  the  diflerence  in  unloading  wave  spi'c-ds  in  l he  three 
1 avers  (2.4/3.  32  3/(i . 4 I /)  . 

A stress  profile  ol  versus  /.-direction  is  shown  in  Fig.  I 1.  Here 
the  existence  ol  lavi'is  is  evident,  the  d i sc<mi  t i nn  i t V ol  radial  stress  at 
(.s'ich  interlaci'  shows  on  1 v app  rox  i nia  I e I v Ihm  .luse  the  inliilace  line  h.is 
lu-en  chosen,  arbitrarily,  to  sit  in  the  upper  rat  In  r I h.in  tin  I owe  i 
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FIG.  lOa  F(7) 


0.0  0.2  a4  ae  o.s  i.o  1.2  1.4 

I ! 1 ~r 1 1 1 1 

PRESSURE 


SOL  10*  Ti^U-DlMt  NSIONAL  OUl-oOlNO  •mvc.  [\  CVLiNDkICAL  COOWOlNATfeS  IrfITH  ft-2  VAftlATlO*^  0*'4LV 

NO  sHfLL*  3 [.AY^■^<^  Triangle  load  <<irH  oepth  vahiatiom 

-•AWflMJ-TEPS  f f On 


0»4  = 

.^UO'JO 

Ul  - 

• j 1 sou 

- 

MAX  PAOIUS  OP  GUIO  = 

7,(<U0 

MAX  OEPTm  of 

GPIO  - ?9.bOO 

)UT=JJT 
NUMncf^ 
NliMHt  H 

c-  V Y 1 (i  T 1 

(;P  POINTS  ON 
OE  P()INIS  ON 

Tl'^t  STEP 
PL  A T c.  In 
SrtELL  IN 

• OP  TO  A MAXIMUM  OF 
P-OlPECIluN  - t> 

2-OlwFCTION  = 30 

aoo  sTtPs 

Ni  |M»iE  P 
1 JMM  E P 

0^'  POINTS  In 
Uh  POIMS  IN 

FIELO  IN 
f III  0 IN 

P-OIPECTION  = 40 

7-niPECTION  = hO 

OUTPUT 

POINTS  ON 

shell  at  1 

I = 

10, 

21f  2di 

29* 

30t  31 

3^-. 

35. 

OUTPUT 

POINTS  IN 

the  fu.lo 

AT  a,j) 

_ 

( 4t  S> ( 

10* 

5) ( 2i « 

51 

{ A*  10)  ( 

lOt 

10) ( 21« 

10) 

p^^SULTS 

COMPi)Tk(>  IN  NuN-DIMh  MS  I ONAL 

FMPM 

LOADING 

ON  INTEPJOW  OF  CAVITY 

“OLLD^  CAVITY.  NO  SHM..L 

L A V 1 T Y 

wAOlJS  = 

1.00 

Cavity 

LLNGIM  = 

1 4 • SO  0 

«A  TFP  I AL 

PftOPFPTIt-S,  FOP 

■»  LAYERS 

L A Y E P 

NLYP  oz 

[ 4 M M n A 

L AM;5uA 

L AmolJA 

NU 

NJ 

Z MAX 

OEnSI TY 

LOAD 

ONLOa J 

■3MEA-J 

LOAD  JULOAO 

i%AllO 

1 

8 .sno 

1 .000 

2.  AOO 

.223 

,467 

3.500 

1.000 

P 

12  -SOO 

1 .OMl 

3.  ^2S 

.'•'JA 

• 2«6 

• 484 

5.500 

1.000 

J 

nO  •‘^U'l 

S,  IaS 

O.  1 A? 

.2d0 

, 364 

29,500 

1 .000 

time  MlSft)»^US  TO  PiuTTPO  EvEWy  /»Th  TIm-:  «;T-3 

V A K I A r>  L f. 
'»-ir  w 

^01 M-  JN  l-ir. 

u 

4 

4 

JOINTS  W T-it  I t i ' 


A 


lOCa  i I :)N 
I J 


1 0 


i U 3 

b 


sOACh  PLJTTr.O  VA*«lAnLl 

*v 

In  InE  rV-OHLCTlOM 
IN  THE  /-OMtCTION 

1 


T I “t  L I NE 

STt  P lUC  a T ion 

150  11 


sigma 

si  uma 

SI 

M 

z 

P 

Z J 1 

J 1 

APPLO 

A 

z 

YZ 

Vf  L 

DISPL 

OliPL 

max 

LOAD 

IT  5 10  T = .l*iO 

FIG.  II 


4i 


FIG.  12 


VELOCITY  TIME  HISTORIES 
CASE  A - NO  SHELL 


FIG.I3  STRESS  PROFILE  IN  THE  Z- DIRECTION 

AT  R=  3.0  , T=  1.5 


CASE  A,  NO  SHELL 


Layer.  The  computation  sclieme  permits  only  one  medium  to  exist  along 
any  grid  line. 


B.  Case  witli  structure  (non-dimensional) 

For  a second  set  of  representative  computations  a structure  was 
added  to  the  material  configuration  of  Case  A.  The  same  loading.  Figs.  lOa, 
10b,  was  applied,  this  time  to  the  interior  of  the  strucuturi-.  Output 
was  obtained  at  the  same  space  points  as  in  Case  A. 

Figure  14  gives  the  initial  POSSl  code  print  out  which  lists  all  the 
parameters,  both  material  and  structural  for  this  case.  The  structural 
parameters  are  the  same  as  those  used  for  the  structure  discussed  in 
Ref.  [4]. 

Figure  15  presents  radial  velocity-time  histories  at  the  same  points 
as  for  Case  A.  The  effect  of  the  m.iss  of  structure,  one  radii  removed 
from  it,  is  apparent  from  the  reversal  of  the  initially  outwardly  directed 

V . 
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C.  Case  with  structure  (dimensional) 


The  input  parameter  chosen  for  this  case  were  similar  to  those  used 
for  the  runs  in  Ref.  [4].  They  are 
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FIG.  15 


VELOCITY  TIME  HISTORIES 
CASE  B - WITH  SHELL 

.8 


r 


At  = 0.01  sec 


Since  the  non-dimensionaiizing  factor  — =1.0  the  velocities  shown  in 
Fig.  15  may  be  interpreted  as  velocity  output  for  this  case.  The  units 
on  the  velocity  scale  in  this  case  would  be  in  inches/sec,  and  the  time 
scale  would  be  in  seconds.  A sample  of  the  computer  output  of  input 


information  is  shown  in  Fig.  16. 
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CONCLUSION 


A description  of  a computer  code  known  as  I’OSSI  has  been  presented. 

The  code  was  written  to  provide  a small,  flexible  tool  for  comput  i nc,  effects 
of  a soil-structure  interaction  problem  in  two-dimensional  (r,  z)  coor- 
dinates, where  tlie  cylindrical  structure  is  buried  in  a type  of  dissipative 
material  wliicli  may  be  layered.  The  POSSI  code  may  be  used  as  a rlieck  code 
for  large  finite  element  or  finite  difference  computations,  or  it  raav  be 
used  to  explore  various  so  i 1 -s true ture  interaction  problems  (within  the 
limitation  of  tlie  type  of  material  presented).  It  may  also  be  used  to 
efficiently  run  CIST  type  problems  for  cases  in  whiclt  the  primary  dissipation 
mechanism  in  tiie  soil  comes  from  pressure-vo 1 ume  hysteresis. 

Results  from  two  sample  problems  have  been  presented.  The  first 
be’'ng  a long  cylindrical  hole,  unlined,  in  a material  of  three  layers.  The 
second  problem  consisted  of  a structurally  lined  long  cylindrical  cavity 
in  tlie  same  three  layered  material.  in  both  cases  results  were  present 
in  non-dimensional  form,  input  and  output  information  from  the  code  was 
provided  for  both.  As  a subcase  to  the  second  the  re.sults  were  also  given 
in  dimensional  form. 

In  any  computations  using  I’OSSI  code  there  is  certain  care  to  tv 
taken.  Boundary  and  corner  points  present  problems  of  approximation  to 
real  situations.  Conclusions  drawn  alnnit  results  .it  such  points  must  be 
more  carefully  checked.  Section  III  of  this  report  lists  some  other  pc'ints 
of  c.ire  which  should  be  taken  when  using  the  i’OSS  I code. 

Results  1 rom  the  two-dimensional  r,  h coordin.ite  problem  for  this 
conf  i gur.i  t ion  h.ive  alreadv  been  presented.  Ref.  |d).  A nvonmietulcd  ste)i 
in  this  work  is  to  comliini’  the  comimter  code  th.it  w.is  written  fi'r 
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Llio  results  of  Ref.  [2]  wiLli  the  I’OSSl  code  to  ])roduce  a t hree-d  iinens  i eiia  1 


code  for  r,  0,  z coordinates.  Sucli  a cod(.'  would  still  be  relativt'ly 
small  and  flexible  enough  to  be  considered  a check  code.  It  would  also 
provide  a code  for  one  of  the  few  three-dimensional  s L rue  t urc‘-med  i urn 
interaction  problems  not  restricted  to  being  in  an  I'lastic  material. 


r 
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Al’PENDiX  A - INPUT  INSTRUCTIONS  I'OR  POSSl  CODE 


All  input  data  ri’quirod  for  tlio  POSSl  code  is  listed  in  Table  6 in 
the  order  in  whicMi  it  is  required.  Variable  names,  deseription  and 
formal  for  eaili  I'ard  set  are  given.  (Variable  names  are  also  defined 
in  Appendix  B.)  Variables  required  for  cboiee  of  options  such  as 
loading  type,  non-dimensional  versus  dimensional  output,  etc.  arc>  also 
1 is  ted. 


Table  7 lists  the  numbering  convent  ion  for  each  variabli’  on  t lu- 
shell  or  in  the  field.  If  printed  output  or  grapbii-  ontpul  is  desiri'd 
tile  (juantity  to  be  oul|)ut  (stress,  velocity  or  displacement)  will  be 
identified  by  this  number  system. 
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TABLE  - DESCRIPTION'  OF  I.WUT  DATA  FOR  the  order  in  which  it  is  required) 
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If  MXCRS  = 0 tifo  following  card  is  to  bu  omitted.  (Tliat  is,  no  time  historic 
of  points  on  the  strurture  are  requested.) 


then  LINE  will  be  a value  fo  the  1 index  in  the  field. 
NRORZ  = 2 (profile  in  the  Z-direction)  then  LINE  will  be 
value  of  the  J index  in  the  field. 


CARD 

SET  VARIABLE  DESCRIPTION  FORMAT 


TABLE 


NU>ffiERING  OF  VARIABLES  FOR  GRAPHIC  OUTPUT  OR  PRINTED  OUTPUT 


Structure 


Field 


VARIABLE 

NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 
9 

1 

2 

3 

4 

5 

6 


VARIABLE 

O Radiul  stress  at  the  structure/material 

interface 

CT^  Longitudinal  Stress 

a Shear  stress 

rz 

V Radial  velocity 
r 

V^  Longitudinal  velocity 

W Radial  displacement 

r 

Longitudinal  displacement 
First  invariant  of  stress 

Oqq  Hoop  stress 

a Radial  stress 

r 

Longitudinal  stress 

a Shear  stress 

rz 

V Radial  velocity 

r ^ 

V^  Longitudinal  velocity 

First  invariant  of  stress 


hi 


LIST  OF  COMMON  VARIABLE  NAMES 


Appendix  B 


AJl 


The  first  stress  invariant,  evaluated  at  points 
in  the  field. 


AJIMX 


Last 


maximum  compressive  value  of 


J 


1 


ALPHP,  ALPHS 
APDSC 


Values  of  u for  plate  and  shell  as  defined  in  Table  1, 
Item  5. 

Load  applied  at  each  point  on  the  Interior  of 
the  shell. 


ARLM1,AR1.P,AZLP 
AZLPl,  AZLS 


Intermediate  quantities  used  in  computation  for 

right  angle  at  shell-plate  interface.  See  Eqs.  (21)-(24) 


CLD 


Coefficient  array  computed  with  loading  P-wave 
parameters. 


CS 

Coef  f icient 

CUN 

Coefficient 

parameters . 

array  computed  with  S-wave  parameters 
array  computed  with  unloading  P-wave 


CNRl CNR12  Coefficients  defined  for  computations  of  motion  of 

shell-plate  structrue  as  follows: 


CNRl  = 

At 

s 

CNR2  = 

8 At 
s 

CNR3  = At 

p 

^ 2 

CNR4  = 

2 

'll  a At 

P P 

r 

CNR5  = 

where  y 

p AR  '"^s 

" D Az  ' 2 

s LC  (1) 

P 

CNR6  = 

LC 

PKTD  7 — 

Y 

CNR8  = B At 
P 

^’'lc-1 

LINK/  — 

1+Y 

CNRll  = 


Ar 


Az 


LC 


CNR12  = 


2[l+(^"')  1 

LC 


DENS 

Density  of  each  layer,  relative  to  a base  value  p . 

DR 

Increment  of  space  in  the  R-direction;  AR. 

DR2.DR3.DR4 

(AR),  (AR),  (AR)^. 

DT 

At  Increment  of  time. 

DZ 

Az  Increment  of  space  in  the  Z-dir-ction.  Possiblv 
different  for  each  layer 

DZ2,DZ3,DZ4 

(AZ)^  (AZ)^,  (AZ)^  . 

f)l 


PRSCSDINO  page  BLaMC-NOT  fllMSD 


l)lRWR,DlR’,vZ,D2RWK 

D2RVJZ,D3RWZ,D4RWZ 

EPL.ESM 

HPL.HSH 

I 

IFLPT 

I LOAD 
ILYR 
IPTS 
IPR 

IR 

IRIZ 
TRSHL 
IRS  PI 
IRZSH 
ISHPR 

ISllPT 

ITGRF 

IZ 

IZSIII, 

IZSPl 

•IFLPT 

• IPR 


Intermediate  quantitie^s  u.sed  in  computing  sliell  velo- 
citie.s  in  the  R and  Z direct  ion.s.  See  Tables  3a, 

3b,  4a,  4b. 

Young's  modulus  for  the  plate  and  shell,  respectively. 

Thickness  of  the  plate  and  shell  respectively. 

Index  in  the  Z-direction  measured  from  tlie  top. 

I-index  of  field  point  at  which  grapii  output  is 
desired . 

Indicator  for  type  of  loading.  See  Table  6,  card  set  2. 
1-inde.x  at  the  bottom  of  eacli  layer. 

Total  number  of  mesh  points  in  the  field,  = IRIZ-IRZSH 

1-index  of  p(jint  in  the  field  at  which  printed 
output  is  desired. 

Maximum  number  of  mesh  points  in  the  R-direction. 

The  product  of  IR  and  IZ. 

Total  number  of  points  on  the  plate. 

= IRSllL+1 

Tlie  product  of  IRSIIL  and  IZSHI. 

Index  of  point  on  shell  at  which  printed  output 
is  desirt'd. 

Total  number  of  points  around  tlu'  cavity  = IRSHI.+  lZSiil.-l  . 
Tinu'  index  increment  for  graplts  of  time  liistories. 

Maximum  number  of  mi>sh  points  in  t lu‘  Z-direction. 

M.iximum  number  of  points  on  tlie  sliell. 

= IZSIlI.+  l 

•1-iiulex  of  field  point  at  wliich  graphic  output  is 
des i red . 

•1-index  of  field  point  at  which  printed  outinit  is 
des i red . 

f)4 


JPTS 


Number  of  mesh  points  in  the  R-clirection  in  ttie 
field  between  the  shell  and  the  outer  boundarv. 
JPTS  = IR-IRSni.. 


JO 

KODES 

l^VMl.D 

LINE 

LR 

LRMX 

LYRSH 

LZ 

LZNLX 

OTPT 

MXGRF 

MXGRS 

MXPRF 

NFLV 

NGRPH 

NI.AYR 

NODIM 

NOR  IT 


Index  in  the  R-direction  measured  from  the  origin. 

Indicator  for  state  of  load,  unload,  reload.  Sec- 
Table  5. 

1 JiiMS  , 1 JiMl'N  Non-dimensional  wave  speeds  for  loading  P-waves, 

unloading  P-waves  and  S-waves,  respectively. 

Index  of  line  along  which  a space  profile  is  to  be 
p 1 o 1 1 ed . 

Array  LR(I)  defines  tlie  number  of  mesh  points  to 
be  included  in  the  integration  in  the  r-direction 
along  each  1-tli  horizontal  line. 

Maximum  values  of  LR. 

Number  of  the  layer  in  which  the-  plate  sits. 

Array  I.Z(.I)  Defines  tfie  numhe-r  oi  mesh  \)oints  to 
be  included  in  the  integral  i oti  in  tlie  z-direction 
along  each  J-tli  vertical  line. 

Maximum  values  of  I.Z . 

Time  stc^p  number  at  whicli  space  profile  is  to  In- 
plotted. 

Total  number  of  time  liistories  desired  for  field  points. 

Total  number  of  time  liistories  desired  for  structure- 
points. 

Total  number  of  space  profiles  desired. 

Field  variable  number  for  wliich  grapiiic  output  is 
desired.  See  Table-  7 . 

Indicator  for  graphic  output  or  not.  Set-  Table  card 
si-t  4. 

Total  mnnber  of  layers. 

Indicator  for  non-dimensional  or  dimensional  results. 

See  Table  b,  card  set  4. 

Indicator  for  right  angle  prt-servat  ion  at  slu-11 
plate-  interface.  See  Tahle  b,  card  set  4. 


(lb 


NOSHL 

NPAIR 

NPRSH 

NPRT 

NRORZ 

NS  HPT 

NSHV 

NT  I ME 
NUED 
NUUN 
NVAR 

OMPL.OMSH 

OM2PL,  OM2SH 

PLDEN.PLNU 

R 

SGR 

SGZ 

SHA.Il 

SHDEN 

SHU  ,SIIE2 

SHNli 


Indicator  for  structure  or  no-structure  lining  tlie 
cavity. 

Total  number  of  field  points  for  which  printed 
output  is  desired. 

Total  number  of  points  around  cavity  at  which 
printed  output  is  desired. 

Time  step  index  increment  at  which  printed  output 
is  desired. 

Indicator  for  direction  in  which  space  profile  is 
to  be  plotted.  See  Table  6,  card  set  11. 

Index  of  tlie  point  on  the  structure  at  which  a 
time  history  is  to  be  plotted. 

Number  of  the  structure  variable  for  which  a time 
history  is  to  be  plotted.  See  Table  . 

Maximum  number  of  time  steps  to  be  computed. 

Poisson  ratio  for  loading  conditions.  See  Table  1, 

Item  4. 

Poisson  ratio  for  unload-reload  conditions.  See  Table  1, 
Item  4. 

Number  of  the  field  variable  for  which  a space 
profile  is  to  be  plotted.  See  Table  7. 

Fundamental  frequencies  for  plate  and  shell 

respect i ve  Iv . ^ ® See  Table  1,  Item  b. 

2 2 
U'  ,li) 

p s 

Density  and  Poisson  ratio  for  plate  material. 

Distance  in  the  r-direction  measured  from  tlie  origin. 
Radial  stress,  at  points  in  the  field. 

Longitudinal  stress,  O^,  at  points  in  the  field. 

First  invariant  at  tlie  structure-field  interface 
Density  of  shell  material. 

Coefficients  for  shell  computations,  SHLl  = l+o^,  SHI. 2-1- 
Poisson  ratio  of  shell  material. 


bb 


SHSGR,SHSGZ,SHSRZ 


measured  in  the  field  at 


Stresses  0 , o , o 

r z rz 


the  St  riieture-material  interface. 

SHVR.SHVZ 

Velocities  of  points  on  the  structure  in  the 
R-  and  Z-d  i rec t ions , respectively. 

SHWR.SHWZ 

Displaci  .nents  of  points  on  the  structure  in  the 
R-  and  Z-d i reit ions , respectively. 

TIME 

Cumulative  time  value  from  start  of  computation 

> 

> 

Velocities  of  points  in  the  field  in  the  R-  and 
Z-direct ions , respectively. 

Z 

Distance  measured  in  the  Z-direction  from  the 
top  boundary. 
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ATTN:  OP  *1S1 

ATTN:  t'P  (APT  !'.>.> !.■ 

ATTN:  Oi’  <»H.\  1 UDK  Smi ! h 

ATTN:  (Y»de  nO-.Ct,  Robert  Pi.i.csi 

ATTN:  Op  qS.'.  1T(  l*ub.!. 

Uhiet  ol  Naval  Kese.ii . h 

N.ivv  Department 


AT  IN:  John  H.iv. 

Oi  k 

AITN:  Nicholas 

Pi*t  roll. 

ATTN:  Ciule 

Th. .na*--  1 

. i.'uiim 

AITN:  Code  Vn.:,. 

Ja.  ol.  ! . 

Warm  i 

ATTN:  Te.  lmli  al 

I ibi at  V 

Ot  J 1*  <*r  in  Ch.ir>'e 
tivil  Ingineeiing  laboi  il-av 
Naval  l onst  Ml.  i ion  I'.atl  ali.'H  (eni.r 
ATTN:  Te.  hni..il  llbi-uv 

ATTN:  R.  J.  ivl.-Mo 

AITN:  SI  m Takah.ishi 
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DK{\\RT>n;NT  OF  THE  NAVY  (ConLin-u-d) 


DEPARTMENT  OF  THE  NA\A  (CotUimu-d) 


Commandani  ot  t!jt*  Marine  Cori  s 
S.ivy  Departr.ent 
ATTN : POM 

CoJiinianJer 

David  W.  Taylor  Naval  Ship  K & I)  Center 
ATTN:  Code  1.42- j,  Library 

Comriand  ing  (leneral 
DfVflopn'.ent  Center 
Fire  Support  Branch 
MCDEC 

ATTN:  CAPT  Hartneady 

AT'TN;  ETC  Oapenski 

Cumm<uider 

Naval  Air  Systems  Command 
Headquarters 

ATTN:  F.  Marquardt 

Commander 

Naval  Electronic  Systeius  Command 
Naval  Electrtmic  Systems  Cmd.  Hqs. 

ATTN:  PMt;  11  7-21 A 

Commanding  Oft'ieer 

Naval  Explosive  Ord.  Disposal  Kac. 

ATTN:  Ci'de  504,  Jim  Petroiiskv 

Commander 

Naval  Eacilit ies  Engineering  Commaiul 
Headi^uarters 

ATTN : lei'hnica  1 Library 

ATTN:  Code  04 B 

ATTN:  Ci)de  OiA 

Comm.inder 

Naval  Ocean  Systems  ilenter 

AI  TN : Tei’hn  ica  1 Li  br.u  > 

Super  int endent  (Code  14241 
Naval  Post  graihia  t e School 

AITN:  tN>de  2124,  Tech.  Kpts.  I.ihrarian 

Directoi 

N.iv.i  1 Re. sea  fell  Labor.!  t i>ry 

AITN;  Code  8*40,  E.  Ri.>senthal 

ATTN:  Code  2h00,  Tech.  Lib. 

Com:Tian'.ler 

N.iv.il  Sea  Systems  Ci>nimand 
Navy  Department 

ATTN:  ORD-0  J j 

ATTN:  Code  0)511 

ATTN:  0Rl>-‘inM.  Lib. 

ATTN:  .SEA-99  nc 

Commamler 

Naval  Ship  Engineering  Center 
DeiMriment  -M  the  N.ivy 

ATTN:  NSEC  h|()5(; 

AITN;  Tei  hni..!!  I Ibr.iry 

Commander 

N.ival  Ship  Rese.ir«h  i>  Developru'nt  Cente? 
Cnd4TW.it  er  Explosive  Rese.ircli  Div  i s l»»n 
AITN:  Te.ljni.  il  I. ibr.iry 


Commander 

Naval  Surface  Weapons  Center 

.vnN:  Code  240,  ('.  I.  Ari»ns4>n 

ATTN:  M.  Kleinerman 

ATT.N:  0.  1..  M.atti-son 

AITN:  Code  W.-\501  , N.ivv  Nuc . Progms.  Oft 

Commander 

Naval  Surface  Weapons  Center 
Da  1 1 1 g re  n La  bo  r a t o ry 

ATTN : Techn ical  1 ibrarv 

Prv'sideni 

N.ival  W.jr  College 

ATTN  : Ti'chni  ca  I I.  ibrarv 

Ciimmande  r 

N.-jval  Weapv'HS  (\*nier 

ATTN:  Code  5j),  IVcii.  I.ib. 

Commanding  ilttic^T 
Nav.il  Weapons  Evaluation  Ea.  i 1 i t v 
ATTN:  Technic.!]  I.ilrarv 

■ATTN:  R.  Hughes 

Di rector 

Strategic  Svstems  Project  Office 
N. i vv  De pa  r t men t 

ATTN;  NSP-272 

.ATTN:  NSP-41.  leeh.  lib. 

AiTN;  NSP-27 I 

.‘>F  JjiA  air  h'oRn: 

Comm.inder 
ADO'M  XPD 

ATTN;  XP 
AITN:  XPOD‘> 

.AF  Ceophvsics  L.ibi*ra t rv , \ES( 

AITN:  Ki  r t‘.  Ih.*mpson 

.\ITN:  Sl'OI  , R.  -a-.n-.-ii  I ib. 

.\i  Institvili-  lO  Techno  1 1'gv , Al 

AITN:  I ibrarv  AKIT.  Hidg.  >.*U.  Ar.  a H 

.Al-  WtMp.m.s  Laboralorv,  \FSC 


A n N : 

DEP. 

I i tm  i . • 1 

. Hr. It  t. Ml 

Al  LN: 

DI  S-S 

, V.  A. 

Pl.irKMJil.'n 

AHN: 

DLS-t 

. R.'bert 

Hennv 

AITN: 

DI  D 

Al  IN: 

Di  S-C 

. Mr.  M. 

l.’i  r 

\rrN: 

sn 

He.ld(|'l.irt  <Tn 

.\ir  Fof.  e Svst4*ms  Comri.iml 

A 1 I N : I'e.  hn  i .‘a  1 I ibi  .it  v 

AHN:  lUiAW 

iriN;  R.  Cross 

C«>mnj.iiJd>T 

\rmanent  Devel  .gneni  io-.t  ..uttj 
AIIN:  h-,  h.  in.iarv 

Al  IN:  ADKKI  - 

Command,  r 
ASD 

AIIN:  I i*.Jm  i . .i  1 1 1 1 i ai  v 


DKPARTMKNT  of  THi;  AIK  FOKa:  (Continued) 


FNKRCV  KKSl.AKCH  4.  ULVl  I Ol'MI  NT 
(Ctnit  imifd  ! 


Assistajjt  Secretary  of  the  Air  Force 


Kesearcli  i Deve U>pnient 

Fn i Vers i t v 

o t C.i  1 11  orn  ia 

Headquarters 

, I’S  Air  Force 

I.awr«.MU  e I. 

i\.’er?:iore  laborat.  rv 

A riN : 

Col  K.  E.  Steere 

ATTN 

D.  T.  Morris,  I.-MO 

ATTN 

Tei-h.  Into.  Dept  . 1 - a 

Deputy  Chiet 

ot  Sta  t f 

ATTN 

Teil  Butkt'vi.  h,  I -J(M> 

Research  Deve lopmeitl 

ATTN 

Larry  W.  Woodruii,  I-'ic 

Headquarters 

, CS  Air  Force 

Al  i N 

R iiTia  rd  G.  Don.' , 1 -qc 

ATTN : 

Col  J.  L.  Gilbert 

Al  I'N 

Jerry  (.KUuir*.Mu 

ATTN 

M.  Fernaiuie;: 

Commander 

ATTN 

J.  K.  Hearst , l.-JDS 

Forei>»n  Technology  Division,  AFSC 

AITN 

Ja<  k Kalin,  L-7 

ATTN: 

KTDP 

ATTN 

Robert  Schock,  1.-4  17 

ATTN : 

I’DBG 

ATTN : 

NICD,  Library 

Sandia  Laboratoj-  ies 

ATiN : 

PDBF,  Mr,  Spring 

Livernon  I 

.aboratorv 

A I TN : 

Hq.  CSAF/IN 

A I'TN : 

IN 

Sand  ia  I aboral or ies 

Al  I N 

Hq.  FSAF/PR 

MTN 

Doc.  Con.  i\»r  . I.  Mi-rritt 

Ai'l  N : 

PKi; 

ATIN 

Iks'.  Con.  fi'r  Luke  J.  \orti\ 

AITN 

Dm-.  Con.  lor  J141,  S.mdia 

Hq.  I'SAF/KD 

AllN 

L.  Hill 

AtTN: 

RDOKM,  Col  S.  C,  Green 

AIT'N 

Doc.  (-'on.  tor  U'.  Koberi'.' 

ATTN : 

KIXIPN,  Maj  F.  Vajda 

AT  I N : 

KDOSM 

FS  Iner,;v  Research  4.  Dev.  Ail'. in. 

AITN: 

RDPS,  Lt  Col  A.  Chiota 

Al  buqtierque 

(iperal  ions  Or  t i.  e 

A ITN : 

RDI’M 

ATTN-: 

Doc.  ton.  :or  .Tsh.  libr.iT'. 

(T>mniander 

FS  Fner -v  R 

esearch  X ]>eV.  Adr;i?j. 

Rome  A i r Devt 

•lopmetu  la-niiT,  AFS( 

Division  i*t 

Headnuar  i ers  S«.  r . i , . 

A7  TN : 

KMiT.D,  Doc.  Library 

Library  Hraiu  Ii  f.-D'.  1 

AilN: 

LMRFC.  U.  N.  Ma i r 

A ! TN : 

Do.-.  Con.  loi  (i.iss  ti.  li.  1 

SAMSO/Dt 

FS  IiH-r/v  Research  i,  lU-v.  .\dmin. 

A 1 TN : 

DLH 

Neyada 

at  ions  i^i  J i.  <. 

ATIN : 

iXic  . Con . for  U.  h . lib. 

SA.MSO,  DY 

AT  r.N : 

D'lS 

Fnion  (arbitie  i.orporat  Ion 

Hoi  it  ielil  N 

at  iona 1 I aboral oi v 

ATTN: 

Doi  . ( .MI . l or  I e«  ii.  lib. 

A 1 I N ; 

MMH 

AITN: 

( i V’  i 1 lU' t . Ri’s . Pfi*  ] . 

A I IN: 

MNNH 

nmi  K i.tn  I.RNMI  Nl  .VT  N(  11  S 

S/VNSn,  RS 

ATTN:  RSS/Col  Donald  Di>wler 

Cent  t a 1 Int ' 

iT  1 i pen.  1 y 

A 1 TN : 

SAMSt)/XK 

NT  I) ’CSI-  NC  .S  bq--. 

A 1 1 N : 

XK I F. 

Depart  i:a-n  t « 

.'1  ( lie  Inl  el  ior 

«.oi‘.:;.ander  in 

( liic  l 

liuteait  O'  T 

lues 

•i » at . i « Air 

(otM'iinil 

AllN: 

1. . h.  1 ib. 

A n N : 

XPf^ 

AllN: 

NR  I -S  I I NT  V),  1 i bra  I \ 

Di  p u t r-  eni  . 

! t hi  Int  1 t j-  J 

It.  'aolo.-,.  . 

1 1 ■'iir  yes 

IMKi.'i  KIMAKlIl  lIKVU.iU'M.St  AIIM 1 M •,  1 KA 1 1 o"; 

Air.: 

.1.  H.  b.  .il'. 

A!  IN: 

( . . i 1 b.  K.i:i  j.-!i 

Divi.i.ai  -’t  -V 

! 1 1 n 1 r*.  Apf>  I i 1 a t i on 

1 .S  Iner/v  R«‘s 

\ \ ‘ \ 

A1  TN: 

U»»  . ( on.  J or  lest  (H  ! i»  e 

A",  . ..  ,1 

b ' eel  1 I 

A .IN; 

Fob.  •:  W.  -la.  ks -n 

I.*»s  Al  ir'ii's  ■>. 

ler.t  i { i.  I.abi'rat  or  v 

AITN: 

( on . tor  . k . :n|i  1 1 I nan 

‘ M * » • 1 ' t Nui  ! 1 a t ‘ 1 1 • r L«  .'ll  1 1 1 1 . n 

/v  1 1 r,  : 

Do-  . t oil.  tor  Kepi  r i lib. 

Nl..  1 • It  Keen  1 1 1 o r •.  ( i >.'>  ion 

A r 1 N : 

I'o»  . » on . tot  A 1 D.i V 1 s 

M ; N : 

Kobi  f 1 F.  in.  '-.in 

DEPARlMtlNT  _p>‘,  CO\TR.\CTORS 


AtTospaCi?  Corporation 

ATTN:  Prom  N.  Matliur 

ATTN:  Larry  Seizor 

2 cy  ATTN:  Tecli.  Info.  Services 

Agbabian  Associates 

ATTN:  M.  Ay.babian 

ATTN:  Carl  Bagge 

Analytic  Services,  Inc. 

ATTN:  Ceorge  Hesselbacher 

Applied  Theory,  Inc. 

2 cy  ATTN:  Tolm  G.  Trul io 

Artec  Associates,  Inc. 

ATTN:  Steven  Gill 

Avco  Research  ^ Systems  Group 
ATTN:  William  Broding 

ATTN:  Research  Lib..  A8)0,  Km.  7201 

Battel le  Memorial  Institute 

ATTN:  Technical  Library 

ATTN:  R.  W.  Klingesmitli 

The  HDM  Corporation 

ATTN:  A.  Lavagnino 

ATTN:  Technical  Library 

The  BDM  Corporation 

ATTN:  Richard  Hensley 

Bel  1 Telephone  Laboratories 
ATTN:  J.  P.  Uliite 

The  Btieing  Company 

ATTN:  R.  H.  Carlson 

ATTN:  Robert  Dyrdahl 

ATTN;  Aerospace  Library 

Brown  Kngineering  Company,  Inc, 

AITN:  Manu  Patel 

California  Institute  of  Technology 
ATT.N:  Thomas  J.  Ahrens 

C.ilifornia  Keseari-h  & Tec!iiu>h>gy,  Inc. 
ATTN:  Technical  Library 

ATTN : Ki’ti  K reyen hagen 

ATTN:  Slieldon  Shuster 

I'niversity  of  Caliloniia 

Berkeley  (iampus  Room  JIH 

Sproiil  Hall 

AITN:  G.  S.u’hman 

(alspan  G*»rporat  ion 

AT'IN:  'lechnii.il  Library 

(enter  tor  Plannlnj;  i.  Ri'Se.ii>!i,  liu  , 

ATTN:  R.  W.  Slmldx-t 

C I ✓ I I / Nim’  I ea  r Sv.st  «‘ms  (!orp. 

AITN:  Robert  (r.iwti»rd 

I’niversitv  ot  Hayttm 

indiisttial  Si  I lit  It  y Super  KL- >0’» 

ATIN:  Hal  lock  C.  Swit i 


DLPAKTMLN'l  tn-  D1  KLNSL  ruNTlC\CTOKS  (Contimied) 

Tniversity  ol  Henver 
C\)  1 o r a do  Sem  i na  r y 
Denver  HeSeareh  Institute 

ATTN:  Sec.  Oflicer  !i>r  1,  Wisi»isV.i 

LG  ii  G,  Inc. 

Albuque’rque  Division 

ATTN:  Technical  I.ibrarv 

Elcctrii-  }\>wer  Research  Institute 
AITN:  George  SI  iter 

Elect  romechanica  I Svs.  of  New  Mexico,  In.., 
ATTN:  R.  A.  'siuink 

Engrg.  Decision  Analysis  (ai..  Inc. 

.ATTN:  R.  P.  Kennedy 

The  Franklin  Institute 

ATTN:  Zenons  Zudans 

Card,  Im-orporated 

ATTN:  G.  L.  Neidliardt 

(k'ne ra  1 Dynamics  (‘orp . 

Pomona  Division 

ATIN:  Keith  Anderson 

iJeneral  Dynamics  Corp. 

Electric  Boat  Division 

ATTN:  Michael  Pakstvs 

Genera  1 Elec  t r ic  Company 

Space  Division 

ValU'v  Forge  Space  Center 

ATTN:  M.  H.  Bonner,  Space  Si  i.  I.ab. 

General  llectric  Companv 

Re-Entry  ik  Environmental  Systems  Div, 

ATTN:  Artluir  L.  Koss 

General  TU'Ctric  (a^mpany 
TEMPi'-t^eni  iT  tor  Ail  vane  ed  St  ud  ies 
ATTN : DAS  I AC 

(a'ni-r.il  Research  (Nirpor.ii  ii»n 

ATTN : Ben  jam  in  A I i-xamler 

Gi'oceiu  er  s I ncorporat  I’d 

ATTN:  Edward  Marram 

H-Tech . La  In' rat  or  ies , Inc  . 

ATTN;  B.  Havtenhaun 

Honevwel  I I tu  orpora  U’ll 
|)el  ensi*  Sv.st  i*ms  Division 
AITN:  T.  N.  Pi-lviv 

I I T Ke.s»-a  r.  h I list  i t ut « 

A TTN : Te.-liiiii  a 1 1 i I'l  ai  v 

ATIN:  K.  I . W.  1.  Ii 

I n i v.  ot  lit  im-i  *.  it  i hi.  ac- 
t'ol  liu’i  of  I nc.  Im  er  Inc 
•ept  . o!  If  el  la  1 s : tl;  1 tu  er  in 
AHN:  Te.l  Hel\f  h}.o 

Inst  i f til , tel  D.  f , II*.,’  \n.i ! •.  ., 

AT  I N : I DA  I I I'l  1 1 I u» , *•'  1 1 , t 
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m PARTMhNT  OK  OKKKNSK  CONTRACTORS  (Cont  itmfd) 

./.  j{.  Wij'gins,  (.'orapany,  Iiu-. 

AI'TN:  Jon  Collins 

Kaiiian  AviUyno 

Division  ot  Kaman  Soiomos  Cot[). 

ATTN : Toolmioal  Library 

AT'I'N:  K.  S.  Crisciono 

ATTN:  Norman  T.  Hobbs 

Kam<in  Soicnoos  C.orporation 

ATTN:  Frank  H.  Shelton 

ATTN:  Paul  A.  Lllis 

A'lTN:  Library 

Kara^o/ian  vS  C.ise 

ATTN:  John  Kara^-o/.  ian 

Loi’klieoJ  Missiles  Spaee  t\>nipanv,  Inc. 

ATTN:  Teehniial  library 

Loeklieeii  Missiles  Spa»,‘e  (Ainipanv,  Ine. 

ATTN:  Ti-m  t.eer:.,  HKIk-  205 

Lovelace  1‘ouiulal  ion  lor  MeJiial  lahu*.  Rsi'li. 
ATTN:  Technical  I ibrarv 

AIIN:  Asst.  Hir.  *>l  Kes.,  Robert  K.  Jones 

Mir  I in  Ma  r ie  1 1 a Aen>sj),iie 
Orlatult)  Hivisii'n 

ATTN:  C.  Kol ieo 

MiTXunu‘1  1 D<*u>;l;is  t!orpi>rat  ion 
All'N:  Robert  W.  Halprin 

MiMillian  Si-ieiu-e  Assov  iates,  Ini’. 

ATTN:  R..berl  Oliver 

Merr  i t I C.ises , 1 ncor  pora  t cil 
ATTN:  J.  L.  Merritt 

A7TN:  'Je<hnical  I.  ibrarv 

Meteorolop.y  Kesearcli,  Inc. 

AITN:  Willi.im  0.  Creen 

The  Mitre  Corporation 
AIIN:  Libr.iry 

Iniversity  ot  New  Mexico 
[)i*pt  . ot  Campus  Security  Police 
A TTN:  (■'.  ii.  I r ianda  ( a I icf  I s 

Natli.m  M.  Newm.irk 

Consul  tiup.  Knpiiieerin^  Services 

111  Ob  A Civil  ITipv  ineer  iny,  Hu  i 1 d in>’. 

Cn iviTs i t V ot  11) ino i s 

ATTN:  Natli.m  M.  Ni-winark 

Pac  i t i ca  Tt-i  hno  1 op.v 
ATTN:  C.  Kent 

AITN:  K.  Iljorlf 

Physi«s  International  Conif*any 


AITN 

Doc  . 

C.m . 

I or 

I .irry 

A.  Bcdirm.inn 

ATTN 

l)c*c  . 

Con . 

1 or 

h.  T. 

M.mre 

Al  IN 

l>o<  . 

(am . 

t ot 

Cove  Vilii  c'tU 

ATIN 

Doc  . 

(ion . 

1 or 

Denn  i s 

nrphal 

Ml  N 

Doc  . 

Con . 

t or 

Rohet  t 

Swill 

AIIN 

Iloi  . 

( on . 

t or 

lech. 

Mh. 

AIIN 

Do,  . 

(ioti . 

t OI 

I red  : 

. Sauet 

A I IN 

!io(  . 

Ct>n . 

t or 

Chat  I. 

H (iod  1 r ev 

OLPAK’IMLNT  0_K  ULFLNSK  CtiNTRACTORS  (Continued) 

Pr«itotype  Deve l<j(>ment  A«s'c»ciales.  Inc. 

ATTN:  '!  . K.  McKinley 

K 6.  D Assoc  i. lies 

AT'I'N:  Henry  Ciu>per 

ATTN:  Harold  L.  Brode 

ATTN:  Arlen  Fields 

ATTN:  Paul  Rausch 

ATTN:  Cyrus  i\  Knowles 

ATTN:  Albert  L.  Latter 

ATTN:  William  B.  Wrigfit,  Jr. 

ATTN:  Robert  Port 

ATTN : Tec  h n i c a 1 Library 

A'f'^N:  Jerry  Carpenter 

AT'lN:  J.  C.  l.ewis 

The  Rand  tiorporation 

ATTN:  Armas  Laupa 

ATTN:  Technical  Library 

ATTN:  C.  C.  Mow 

Scieni'e  Appl  ications,  Inc. 

ATTN:  I).  F.  Maxwell 

AT'I'N:  David  Bernstein 

Science  Applications,  Inc. 

ATTN : Tec  hn  i ca  1 1.  i bra  ry 

Sc-  ience  Appl  icat  ions,  I tu- . 

AT'I'N:  James  Cramc-v 

Science  Applications,  Inc. 

A'l'TN:  Steve  Oston 

Science  A[>pl  icat  ions,  Inc'orpor.it  ed 
ATTN:  William  M.  L.ivson 

AT'I'N:  Burt  Chambers 

Sc'ut  hwest  Researc  h 1 nst  i t ul  e 
ATTN;  A.  H.  Wen/el 
ATTN:  Wilt  red  i.  Baker 

Stanford  Rosearch  Institute 

ATTN:  Cec'fp.c'  K.  Abt.ihamsou 

Systc-ms,  Scieiu  t*  6.  Softw.iie,  liu  . 

ATTN:  Kt'hei  t Sedpewii  k 

ATTN:  llu'mas  D.  Rim-v 

ATTN:  Dc>nald  K.  (irine 

ATTN:  Tc’d  Chc-rrs’ 

AT  TN  : I i-c  till  i . .1 1 I 1 hi  .irv 

Terr. I Tek  , I in  . 

A TTN  : 1 ec  Imic  al  1 i hi  ai  v 

AITN:  Siduev  t'leen 

AITN:  A.  11.  t.mes 

Tet  ra  !'i-c  h . , 1 lu  . 

ATTN:  I 1-S.in  Hwanp 

AT  I'N : T i‘«  fm  i i .i  I I i hi  a i v 

Tex. IS  A .S  M I'niversitv  System 
lex, IS  A N M Rese.iich  louiulalion 
ATIN:  Hariv  Covie 

1 RW  Delense  N :-p.l« f Svs.  l-li'ilp 

ATTN:  lech.  Into,  t mt  ei  1“  H» 

ATTN:  Nc'im  1 ipnei 

. N AITN:  P.-f  ot  K.  Dai.  K1  '.’1  ;0 

AT'IN:  Donal.t  lotto. a.  Kl-.’l.. 

AIIN:  i’l  av  In  Bhut  t a.  R1  1 Hi  • 


r 


i 

Dt P AR'I’MKN'r  0F_  DEFENSE  CONTRACTORS  ( CoiU  i niit-d ) 

TRW  Defense  Space  Sys.  Croup 
San  Bernardino  Operations 

ATTN;  F.  Y.  Wong,  527/712 

TRW  Defense  & Space  Sys.  Group 
ATTN:  Gregory  D.  HuLcher 

^ Universal  Analytics,  Inc. 

i ATTN:  F.  I.  Field 

i 

J URS  Research  Company 

^ ATTN:  Technical  Library 

, ATTN:  Ruth  Schneider 

The  Eric  H.  Wang 

Civil  Engineering  Research  Fac. 

University  Station 

ATTN:  Larry  Bickle 

ATTN:  Neal  Baum 

I 


DEPARTMENT 

OF  DLFiiNSK  CONTRACTORS  (Continued)  \ 

Washington 

State  University 

Administ  rat ive  Organizat ion 

ATTN 

Arthur  Miles  lU>Iu>rf  for  George  Dtiv.il 

We  id  1 inger 

Assoc.  Consult ing 

Engineers 

ATTN 

J.  M.  McCormick 

ATTN 

Melvin  L.  Baron 

ATTN 

A.  T.  Matthews 

ATTN 

H.  H.  Bleich 

Weidl  inger 

Assoc.  Consulting 

Engineers 

ATTN 

J.  Iscnberg 

West  ingliouse  Electric  Corp. 
Marine  Division 

ATTN 

W.  A.  Volz 

; 

1 


1 


